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灾难 性 测 光 红 移 对 重子 声波 振荡 和 弱 引 力 透 镜 
限制 暗 能 量 状态 方程 参数 的 影响 ， 
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摘要 : 研究 了 存在 灾难 性 测 光 红 移 误差 时 重子 声波 振荡 (Baryon Acoustic Oscillations, BAO) 
和 弱 引 力 透 镜 ( Weak Lensing, WL) 的 互补 效应 ， 以 及 灾难 性 测 光 红 移 误差 对 限制 暗 能 量 状 
态 方程 参数 的 影响 。 针 对 类 大 口径 全 天 巡视 望远镜 (Large Synoptic Survey Telescope, LSST) 的 
RRRA, EI T E zzp EL AUL) AE TF f (BR) 局 域 性 分 布 的 灾难 性 测 光 红 移 误差 
模型 ， 并 分 别针 对 左上 角 和 右 下 角 的 灾难 性 测 光 红 移 ， 利 用 费 舍 尔 矩 阵 (Fisher Matrix) 分别 
估计 它们 对 重子 声波 振荡 、 弱 引力 透镜 和 联合 重子 声波 振荡 + 弱 引 力 透 镜 (BAO+WL) 限 制 暗 
能 量 状态 方程 参数 的 影响 。 若 拟 合 模型 没有 包括 实际 存在 的 灾难 性 测 光 红 移 误差 ， 则 左上 角 和 
右 下 角 的 灾难 性 测 光 红 移 造成 的 偏差 并 非 总 是 同 号 的 。 重 子 声 波 振 荡 受 灾难 性 测 光 红 移 误差 的 影 
响 最 小 ， 对 于 总 比例 为 下 =0. 02 的 灾难 性 测 光 红 移 ， 最 大 的 影响 也 只 是 左上 和 角 和 右 下 角 两 部 分 对 
W0 造 成 的 偏差 ， 相 对 偏差 约 为 30%。 但 在 对 弱 引 力 透 镜 和 联合 重子 声波 振荡 + 弱 引 力 透 镜 的 影响 
上 ， 通 常 都 有 几 倍 的 偏差 ， 因 此 灾难 性 测 光 红 移 误差 的 影响 不 可 忽略 。 另 外 ， 关 于 弱 引 力 透镜 对 
wo 的 限制 ， 左 上 角 和 右 下 角 造 成 的 偏差 大 小 几乎 相等 ， 但 符号 相反 ， 从 而 导致 总 体 影响 反而 变 得 
很 小 。 当 拟 合 模 型 包括 灾难 性 测 光 红 移 误差 在 内 后 ， 虽 然 测 光 红 移 误 差分 布 模型 多 了 45 个 自由 
度 ， 但 在 同样 大 小 的 训练 集 的 情况 下 ， 重 子 声波 振荡 和 弱 引 力 透 镜 的 互补 效应 仍然 很 强 。 在 此 条 
件 下 ， 瞳 能 量 状态 方程 参数 的 误差 并 没有 太 大 的 增加 。 特 别 是 重子 声波 振荡 的 限制 结果 增加 
量 少 于 1%。 弱 引力 透镜 对 wo 和 w, 的 限制 误差 分 别 增加 了 14% 左右 (UL+BR) 和 6% 左 右 (UL+ 
BR) ， 而 联合 重子 声波 振荡 + 弱 引 力 透 镜 对 w Few, 的 限制 误差 都 只 增加 5% 左 右 (UL+BR) 。 
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弱 引 力 透 镜 被 认为 是 探测 暗 能 量 的 一 种 强 有 力 的 工具 ， 它 通过 探测 遥远 星系 图 像 的 切 变 来 研究 暗 
能 量 的 性 质 及 其 他 宇宙 学 参数 。 很 多 大 型 巡天 望远镜 (如 LSST? 、EUCLID2 、WFIRSTS ) 把 弱 引 力 透 
镜 作为 探测 暗 能 量 的 重要 手段 。 由 于 弱 引 力 透 镜 的 信号 极其 微弱 ， 为 了 提取 有 用 信息 ， 未 来 的 大 型 弱 
引力 透镜 项 目 需要 探测 大 量 的 星系 红 移 样本 ， 以 至 于 探测 全 部 样本 的 光谱 红 移 在 时 间 上 是 不 可 能 的 。 
因此 ， 必 须 用 多 色 测 光 的 方法 估计 红 移 ， 而 测 光 红 移 通过 多 波段 成 像 ， 利 用 经 验 公 式 02 、 能 谱 分 布 
(Spectral Energy Distribution ，SED ) BARE! 、 人 工 神经 网 络 趾 等 方法 估计 红 移 。 对 宇宙 学 有 用 的 是 
星系 的 红 移 分 布 ， 而 不 是 每 个 星系 的 红 移 。 通 过 抽样 测量 星系 的 真实 红 移 或 者 利用 宇宙 学 方法 可 以 校正 
测 光 红 移 误差 分 布 :5 5 。 但 是 红 移 分 布 的 不 确定 性 不 可 避免 地 传递 到 暗 能 量 状态 方程 参数 及 其 他 宇宙 学 
BB EO! 。 研 究 表明 ， 若 所 有 在 真实 红 移 区 间 为 0. 1 的 星系 的 测 光 红 移 误差 分 布 都 可 用 一 个 平均 值 在 
真实 红 移 附近 的 高 斯 分 布 近似 ， 则 需要 10° 到 10° 的 光谱 红 移 校正 ， 才 使 得 暗 能 量 的 误差 放大 限制 在 1.5 
ARDY?! 。 其 次 ， 由 于 星系 能 谱 分 布 中 的 谱 线 特征 信息 混淆 或 者 波段 缺乏 特征 信息 ， 导 致 估计 的 测 光 
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红 移 对 真实 红 移 出 现 较 大 的 偏离 Se， 即 所 谓 灾难 性 测 光 红 移 。 特 别 是 红 移 约 小 于 1 的 星系 能 谱 分 布 中 
的 巴尔 英 跳 变 ( Balmer break) 和 红 移 约 大 于 2 的 星系 能 谱 分 布 中 的 拉 曼 跳 变 ( Lyman break) 的 混淆 导致 
局 域 性 的 灾难 性 测 光 红 移 吕 ] 。 灾 难 性 测 光 红 移 的 存在 使 弱 引 力 透 镜 对 暗 能 量 状态 方程 参数 的 估计 产 
A dP! 。 对 于 灾难 性 的 测 光 红 移 ， 可 以 通过 剔除 一 些 灾 难 性 红 移 可 能 性 较 大 的 样本 ， 从 而 在 一 
定 程度 上 减少 灾难 性 红 移 导 致 的 系统 误差 的 同时 又 不 会 过 多 地 增加 因 减 少 样本 而 导致 的 统计 误差 。 
星系 功率 谱 (或 星系 关联 函数 ) 中 的 重子 声波 振荡 包含 了 在 最 后 散射 截面 的 声波 视界 (标准 尺 ) 的 
信息 。 利 用 标准 尺 可 以 测量 哈 勃 参数 - 红 移 关系 及 角 直 径 距 离 - 红 移 的 关系 '5] 。 在 多 色 成 像 巡 天 中 ， 
仅 靠 重子 声波 振荡 测量 的 多 个 红 移 处 的 角 直 径 距 离 对 暗 能 量 状态 方程 参数 的 限制 有 较 大 的 简 并 性 59 。 
此 外 ， 由 于 星系 偏 祖 因 子 的 不 确定 性 ， 使 得 难以 从 星系 角 功 率 谱 中 抽取 更 多 的 信息 :58 。 
由 于 重子 声波 振荡 和 弱 引 力 透 镜 的 互补 效应 ， 当 把 星系 与 弱 引 力 透 镜 的 交叉 功率 谱 也 考虑 进来 
时 ， 在 没有 灾难 性 测 光 红 移 的 情况 下 ， 一 方面 弱 引 力 透 镜 可 以 帮助 校正 星系 偏 祖 因 子 ， 另 一 方面 不 同 的 
测 光 红 移 区 间 的 星系 交叉 角 功 率 谱 对 测 光 红 移 的 误差 分 布 非常 敏感 ， 从 而 对 测 光 红 移 误差 分 布 参数 具 
有 很 强 的 自 校 正 作 用 ， 联 合 重子 声波 振荡 和 弱 引 力 透镜 也 能 使 测 光 红 移 偏差 gz 得 到 更 强 的 限制 '” 。 
本 文 研究 由 于 存在 局 域 性 的 灾难 性 测 光 红 移 但 没有 把 它们 纳入 拟 合 的 模型 而 导致 暗 能 量 状态 方程 
参数 在 重子 声波 振荡 和 弱 引 力 透 镜 单独 以 及 联合 限制 下 的 估计 偏差 。 然 后 探讨 把 灾难 性 测 光 红 移 纳入 
> 拟 合 模型 后 重子 声波 振荡 和 弱 引 力 透 镜 的 互补 效应 和 对 暗 能 量 状态 方程 参数 的 限制 能 
0 本 文 第 1 节 主 要 介绍 星系 功率 谱 和 和 弱 引力 透镜 功率 谱 ， 包 括 局 域 竹 的 灾难 性 测 光 红 移 在 内 的 测 光 红 移 误 
= 差分 布 的 模型 ， 以 及 误差 估计 的 方法 一 一 费 舍 尔 矩阵 ; 第 2 节 给 出 本 研究 的 预测 结果 ; 第 3 节 是 本 文 的 绪论。 


O 1 方法 


一 1.1 观测 量 ;功率 谱 
以 星系 的 角 功 率 谱 、 弱 引力 透镜 的 角 功率 谱 以 及 它们 的 交叉 功率 谱 作为 观测 量 。 在 林 伯 (Limber) 
近似 -20 下 ， 星 系 和 弱 引 力 透镜 的 角 功 率 谱 可 以 写 为 中 


PEO =T aO D(a) NOWO ACL, 2) , (1) 


HP, cE; HC), D(a) M ALC, z) 分 别 是 红 移 为 处 的 哈 蔓 参 数 、 共 动 角 直 径 距离 和 无 量 岗 
的 物质 功率 谱 ， 窗 函 数 W! (2) 可 写 为 


= n, (2) 
二 b(z) 一 一 X=g 
Q Fia n, 5 
if l 3 Q, H D, (z) k ,ni(z") D, (z, z') _ ( ) 
2 H(z) ac Z 元 Dy’) Á 


其 中 ，g My 分 别 代表 重子 声波 振荡 和 弱 引 力 透 镜 ; b(z) 是 线性 的 星系 成 团 性 偏 祖 因 子 ; 2, 和 H, 4 
别 是 红 移 为 0 处 的 物质 成 分 比 和 哈 勃 常数 ;wu 是 宇宙 尺度 因子 ; n,(z) 和 元 分 别 是 第 ;个 测 光 红 移 区 间 
的 真实 红 移 分 布 和 星系 数目 。 重 子 声波 振荡 和 弱 引 力 透 镜 的 分 区 办 法 可 以 相互 独立 。 
在 实际 观测 中 ， 由 于 星系 角 功 率 谱 受到 星系 散 粒 (形状 ) 噪 声 元 ! 的 影响 ， 弱 引力 透镜 角 功 率 谱 则 
受到 固有 形状 噪声 y2, 元 1 的 影响 ， 因 而 观测 的 功率 谱 为 
BY (1) = po (1) + aE, (3) 


其 中 ， 8k All 2 是 克 罗 内 克 (Kronecker delta ) 函数 。 对 于 星系 角 功 率 谱 ， Lim =l; 对 于 弱 引 力 透 镑 ，y，~0.2。 
本 文 用 CAMBI 2 输出 不 同 红 移 处 的 非 线性 物质 功率 谱 (基于 Halofit 模型 [ 当 的 修正 "1 )， 具体 
使 用 参考 CAMB 主页 。 
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1.2 模型 与 巡天 属性 

模型 参数 包括 10 个 宇宙 学 参数 ，31 个 星系 偏 祖 因 子 ，60 个 测 光 红 移 误差 分 布 参数 。 当 把 灾难 性 
测 光 红 移 加 入 拟 合 模型 时 ， 会 增加 45 个 灾难 性 测 光 红 移 误差 分 布 参数 。 

1.2.1 模型 

在 没有 灾难 性 测 光 红 移 的 情况 下 ， 每 一 个 测 光 红 移 区 间 的 真实 红 移 分 布 用 高 斯 函数 描述 。 如 果 存 
在 灾难 性 测 光 红 移 ， 则 可 相应 多 加 一 个 高 斯 函数 描述 灾难 性 测 光 红 移 的 真实 红 移 分 布 。 

假设 巡天 的 测 光 红 移 分 布 是 n(z)， 对 测 光 红 移 空间 平均 分 成 N 个 区 间 ， 在 没有 灾难 性 测 光 红 
移 的 情况 下 ， 第 i 个 区 间 星 系 的 真实 红 移 分 布 可 用 高 斯 分 布 近 似 : 

1 | (2 = ee | 
-一 一 一 exzp| -一 |， 
VITO core. ; 20 core, i 
其 中 ，z,,,;、o60 ;分 别 是 第 i 个 区 间 的 平均 红 移 以 及 该 分 布 的 标准 差 。 

存在 灾难 性 测 光 红 移 的 情况 下 ， 任 意 区 间 的 真实 红 移 分 布 可 用 两 个 高 斯 分 布 的 又 加 近似 ， 分 别 用 
于 描述 核心 (真实 红 移 附近 ) 部 分 和 灾难 性 测 光 红 移 部 分 的 真实 红 移 分 布 : 


(4) 


Noore, i (z) > 


Not, i(z) =(1 -=f;) Nore, i(z) +My. i(z) ’ (5) 
其 中 , f E i PA DOE PEM CAA% H AK A A ECB; (2) 是 该 区 间 的 灾难 性 测 光 
红 移 的 真实 红 移 分 布 : 
(z — Zeat D 
New, i(2) = exp - | ; (6) 
2T Oon, i 20 cat, i 


HP, zai Ooa ;分别 是 该 分 布 的 平均 值 和 标准 差 。 
参照 [9, 19] ， 把 测 光 红 移 为 0 到 3 的 区 间 均 匀 分 成 30 个 子 区 间 ， 即 N=30， 每 个 子 区 间 宽 为 0. 1。 


zai 和 0; 的 基准 值 取 : ; 
Zox. 1 = 0. li + 0.05, (7) a. 
Ooore, i 50.0501 + z0, i) - (8) 2.5 Po 

i 值 取 0 到 29， 即 每 个 区 间 的 中 心 值 作 为 该 区 间 

的 星系 真实 红 移 分 布 的 平均 值 。 把 这 30 个 区 间 a 


作为 重子 声波 振荡 的 区 间 (0-29)， 然 后 每 5 个 
区 间 合 并 成 一 个 区 间作 为 弱 引 力 透镜 区 间 ， 即 共 me LS Ee 
有 6 个 弱 引 力 透 镜 区 间 (0-5) 。 弱 引力 透镜 区 间 
的 星系 分 布 和 数目 可 以 通过 组 成 它们 的 重子 声波 lr 
振荡 区 间 的 星系 分 布 和 数目 得 到 。 描 述 核 心 部 分 vs 
的 真实 红 移 分 布 的 参数 为 30x2=60 个 。 ss Oe 
由 于 光谱 特征 的 混淆 或 者 光谱 缺乏 信息 量 等 “3 
原因 ， 测 光 红 移 的 估计 会 出 现 灾难 性 的 错误 ， 特 D 
别 是 当 光谱 上 的 拉 曼 跳 变 和 巴尔 莫 跳 变 混淆 时 ， z 
会 在 z-z, 面 的 左上 角 和 右 下 角 分 别 出 现 局 域 性 。 ”图 1 灾难 性 测 光 红 移 误差 分 布 示意 图 。 测 光 红 移 从 0 
分 布 的 灾难 性 测 光 红 移 中 (下 文 将 用 UL 和 BR 到 3 分 成 30 个 区 间 ， 灾 难 性 测 光 红 移 分 布 在 右 下 
分 别 指 代 左 上 角 和 右 下 角 出 现 的 灾难 性 测 光 红 GL ae sd. 
y J a 个 区 间 ) EW ADEE eH RAEN T A 
移 ) 。 为 构造 在 =-zw 面 的 左上 角 和 右 下 角 分 别 出 了 
现 局 域 性 分 布 的 灾难 性 测 光 红 移 ， 令 在 前 5 个 区 Fig.1 Schematic diagram of the catastrophic photo-z error 
间 和 最 后 10 个 区 间 才 存在 灾难 性 测 光 红 移 ， 而 


distribution. The Photo-z space is divided into 30 bins 


第 5 到 第 19 个 区 间 没 有 灾难 性 测 光 红 移 ， 见 图 from 0 to 3. The catastrophic photo-z is distributed at 
1, 前 5 个 区 间 的 参数 基准 值 统一 为 Zari 2.5, the bottom-right Co" 三 和 bin) and the upper-left 


CE 0.5(i=0 到 4), 最 后 10 个 区 间 的 参数 基 (20"-29" bin) of the z-z,, plane 
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hinay h£ 
\ HT IGATV 


ERAT 


4 Kx 研究 与 R 13 卷 


准 值 统一 为 z,, ,=0.4，o, ;=0.3(i=20 到 29)。 对 所 有 存在 灾难 性 测 光 红 移 的 区 间 / 取 相 同 值 ， 并 
根据 给 定 的 灾难 性 测 光 红 移 的 总 比值 到 确定。 用 N; 表示 测 光 红 移 分 布 归 一 化 为 1 时 任 一 测 光 红 
移 区 间 的 星系 数目 ， DN, 表示 包含 灾难 性 测 光 红 移 的 区 间 的 星系 数目 总 和 ， 则 f 和 F, 的 关系 满足 
F, =f} N, 。 本 文采 用 (13) 式 作为 星系 的 测 光 红 移 分 布 ， 经 计算 /和 F, 的 关系 满足 /=3.65F,。 描 
述 灾 难 性 测 光 红 移 的 参数 共有 15 x3=45 个 。 

参照 文 [ 19] 的 方法 ， 在 真实 红 移 为 0 到 3 分 成 30 个 区 间 ， 每 个 节点 对 应 一 个 星系 偏 祖 因子 参数 ， 
共 31 个 ， 其 他 红 移 上 的 偏 祖 因 子 通过 线性 插值 得 到 。 星 系 偏 祖 因 子 的 基准 值 随 真 实 红 移 的 变化 关系 
为 6=1+0. 84z0 , 

假设 每 个 重子 声波 振荡 区 间 用 mu 个 光谱 红 移 校正 。 测 光 红 移 模型 中 参数 的 先 验 误差 与 ,的 关系 为 '™ 


Azero, i = Cem, i/ Vm X (1 =F) ， (9) 
Agi. i E AZo, Wl, (10) 

Azen, 1 = Fon, / Vie Xf » (11) 
Ac ;= Aza of 2 ， (12) 


本 文 取 n,=400，Alnf=100%。 星 系 偏 祖 因 子 先 验 误差 取 Alnb=15%。 
宇宙 学 参数 取 10 个 ， 它 们 分 别 是 暗 能 量 状态 方程 参数 w Mu, ETHE w,， 有 上 暗物质 密度 w.， 
曲率 的 等 价 成 分 比 Q,.， 约 化 的 哈 勃 常数 刀 ， 原 初 氨 丰 度 Y， 谱 指数 n.， 谱 指数 的 移动 (running)al 定 
义 同 [19] )， 原初 曲 率 功率 谱 在 =0.05Mp c 处 的 幅 值 Ai 。z 和 ww, 是 描述 随 红 移 演化 的 暗 能 量 状态 
WHE: w=wtw, (1-a) 。 这 些 参数 基准 值 和 先 验 误差 主要 参考 了 PLANK2013 年 的 结果 '”i ， 具 体 见 表 1。 
表 1 宇宙 学 参数 基准 值 及 其 先 验 误差 


Table 1 Fiducial cosmological parameters and their prior errors 


Wo w, w, w, 2, h Y, n, a Ai 
基准 值 -1 0 0. 022 0. 12 0 0. 673 0. 248 0.96 -0.013 2. 196E-9 
先 验 误差 0.5 2.5 2.8E-4 2.7E-3 0.023 0.012 1.2F-4 7.3E-3 9.0E-3 5. 6E-11 


1.2.2 巡天 属性 
本 文 依据 大 口径 全 天 巡视 望远镜 的 巡天 能 力 设 定 巡 天 参数 。 取 星系 的 测 光 红 移 分 布 为 
Ny, & z exp( = Zn/ 20) (13) 


取 z=0.5， 巡 天 面积 取 area=20 000 deg’, HAM MBER N, = 50 arcmin’, 
参照 [ 28 ] ， 本 文 取 


其 中 ，z, 为 弱 引 力 透 镜 区 间 中 心 值 。 
1.3. R&R (Fisher Matrix) 

用 N 和 NN, 分 别 表示 重子 声波 振荡 区 间 的 数目 和 弱 引 力 透 镜 区 间 的 数目 ， 用 ,表示 它们 的 数目 
ZA, HN, =N+N =36。 

参数 po 的 费 舍 尔 和 矩阵 可 表示 为 2 

a07(1) ,6000) 
> Po Pg 
其 中 ， 0(1) 是 一 列 由 独立 的 观测 量 PEC) 构成 的 列 向 量 ， 共 有 NN, + 1)7Z2 个 元 素 。 重 子 声波 振荡 ， 
弱 引 力 透镜 和 联合 重子 声波 振荡 + 弱 引 力 透 镜 分 别 对 应 于 N = N,, N, 和 N,,， 即 

Oardan() =P,(D) (N>m=n=0), (16) 

其 中 ，m(n) 是 重子 声波 振荡 区 间或 (和 ) 弱 引 力 透 镜 区 间 ( 一 起 ) 排列 的 序列 号 ; P,,, 是 相应 的 功率 谱 户 " 。 
m 和 满足 关系 : 


Yms = 0. 18 + 0. 042z; , (14) 


Fig = 


+ das (15) 
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7 i if X— BAO 
| GN, if XOWL ’ c1) 
C 是 观测 量 0 的 协 方差 矩阵 ; C7 是 C RAEE, EIERE AEF, C 的 元 素 是 


[B, BP, + Py, Pal, (18) 


Cw = (2L + 1) As 
E, u=m(m+1)/2+n; v=r(r+1)/2+s; Py EMBRAER ate LEB 
ERI JT Fe FEE SAE) Be a OR EE aa, AS CTT EE BR il ee 2 A 
缘 化 后 的 误差 ， 即 o(p,) = [Fee o 
假设 观测 量 存在 系统 误差 ， 则 拟 合 的 参数 会 出 现 偏离 ， 可 以 用 费 舍 尔 矩 阵 估计 系统 误差 给 参数 带 
来 的 偏差 |8p。| 0, 


00, 
Op, = [Fle > AO, LG “J — ? 
1 pg 


其 中 ，AO, 是 由 于 系统 误差 给 观测 量 0, 带 来 的 偏差 ， 下 标 遵 循 爱 因 斯 坦 求 和 约定 。 


2 预 测 


为 避免 可 能 存在 的 暗 能 量 成 团 性 效应 ， 本 文 忽略 大 尺度 的 星系 及 其 图 像 切 变 统计 信息 ， 同 时 为 避 
免 小 尺度 非 线性 效应 和 重子 效应 “i ， 对 于 弱 引 力 透 镜 取 40<!<2000， 对 于 重子 声波 振荡 取 40<1<3000。 
且 重 子 声波 振荡 还 要 求 功率 谱 AS (k; z) <0.4， 在 实际 处 理 时 z 取 相 应 重子 声波 振荡 区 间 的 中 心 值 。 
2.1 由 灾难 性 测 光 红 移 导致 的 暗 能 量 状态 方程 参数 的 系统 偏差 

K 2 是 没有 灾难 性 测 光 红 移 误差 时 对 暗 能 量 状态 方程 参数 的 限制 情况 以 及 由 F, = 0. 02 的 灾难 性 
测 光 红 移 误差 导致 的 暗 能 量 状 态 方程 参数 的 偏差 。 其 中 ，o (p) 是 没有 灾难 性 测 光 红 移 误差 时 的 限制 
情况 ; 5, 是 灾难 性 测 光 红 移 导致 的 偏差 ; 6,/o(p) 即 两 者 比值 。6, 和 6,/o(p) 下 分 别 有 3 列 ， 分 别 对 
应 于 图 1 中 左上 角 出 现 的 灾难 性 测 光 红 移 ， 右 下 角 出 现 的 灾难 性 测 光 红 移 和 左上 角 + 右 下 角 出 现 的 灾 
难 性 测 光 红 移 的 情况 。 注 意 F 是 指 包括 左上 角 出 现 的 灾难 性 测 光 红 移 和 右 下 角 出 现 的 灾难 性 测 光 红 
移 在 内 的 总 比值 ， 当 只 有 左上 和 角 出 现 的 灾难 性 测 光 红 移 或 右 下 角 出 现 的 灾难 性 测 光 红 移 误差 时 ， 相 应 
的 f 仍然 满足 f=3. O5F,. 

从 表 2 和 图 2 可 以 看 出 ， 左 上 角 出 现 的 灾难 性 测 光 红 移 和 右 下 角 出 现 的 灾难 性 测 光 红 移 的 影响 有 
些 是 正 负 反 号 的 (重子 声波 振荡 对 w,， 弱 引力 透镜 对 w,。， 重 子 声波 振荡 + 弱 引 力 透 镜 对 wo ) ， 特 别 在 
弱 引 力 透 镜 限制 w 的 影响 上 ， 左 上 角 出 现 的 灾难 性 测 光 红 移 和 右 下 角 出 现 的 灾难 性 测 光 红 移 造成 的 
偏差 大 小 几乎 相等 ， 但 符号 却 相 反 ， 从 而 导致 总 体 的 影响 反而 变 得 很 小 。 重 子 声波 振荡 受灾 难 性 测 光 
红 移 误差 的 影响 最 小 ( 见 表 2 和 图 2)， 对 于 f=0.02， 最 大 的 影响 是 左上 角 和 右 下 角 出 现 的 灾难 性 测 
光 红 移 两 部 分 对 wo 的 偏差 (相对 偏差 为 31%， 见 表 2)。 

表 2 了 暗 能 量 状态 方程 参数 因 灾 难 性 测 光 红 移 ( ,=0. 02) 导致 的 系统 偏差 
Table 2 Biases on the dark energy EOS parameters caused by the catastrophic photo-z error (F,=0. 02) 


dp ôp/o (p) 
o(p) 

UL’ BR? BOTH“ UL BR BOTH 

Bie Wo 0. 1188 -0. 031 -0. 0058 -0. 037 -0. 26 -0. 049 -0.31 
w, 0. 4191 -0. 0097 +0. 04 +0. 030 -0. 023 +0. 095 +0. 072 
WL Wo 0. 088 -0.75 +0. 69 一 0. 058 -8.5 +7. 8 —0. 66 
w, 0. 340 -1.1 -0. 42 -1.4 -3.1 -1.2 -4.2 

Wo 0. 064 -0.31 +11 +0. 75 -4.8 +17 +12 

BAO+WL 
0. 2081 -1.0 -2.7 -3.6 -4.9 -13 -17 


a) include UL only; b) include BR only; c) include UL+BR 
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图 2 6(wo)/o(wo) 以 及 6(w,)/o(w,) 随 下 的 交 化。 左边 3 幅 图 从 上 到 下 分 别 是 图 1 中 左上 角 、 右 下 角 以 及 这 两 部 分 
灾难 性 测 光 红 移 误差 对 w 造成 的 相对 偏差 。 右 边 3 幅 图 同 左 ， 但 是 对 zw, 造成 的 相对 偏差 

Fig.2 The variations of 6(w,)/o(w,) and 6(w,)/o(w,) with F,. The three panels on the left are the impact of the upper-left 

(UL), bottom-right (BR) and both the two parts (UL+BR) catastrophic photo-z in Fig. 1 on wọ. The right panels are 


the same as the left but for impact on w, 


DEVE WM GEL BS wre Fe OY E TE IRA AXI | TA BEN A JET ae h Fal a eB aT K 
DUCE THB aN RGSS , TCM TEL he Boe AY DT) EE, PA S| 7 a A PL Be 
荡 受 灾难 性 测 光 红 移 误差 的 影响 更 明显 ， 从 而 直接 导致 弱 引 力 透 镜 的 功率 谱 比重 子 声波 振荡 受灾 难 性 
测 光 红 移 误差 的 影响 更 明显 。 

2.2 把 灾难 性 测 光 红 移 的 误差 分 布 加 入 拟 合 模型 中 时 暗 能 量 状态 方程 参数 的 限制 情况 

在 没有 灾难 性 测 光 红 移 误差 的 情况 下 ， 由 于 弱 引 力 透 镜 不 受 星 系 偏 祖 因子 的 影响 ， 联 合 重子 声波 
振荡 可 以 有 效 地 提高 对 星系 偏 祖 因 子 的 限制 能 力 '” 。 让 测 光 红 移 误差 模型 包括 灾难 性 测 光 红 移 误差 ， 
结果 发 现 ， 对 于 =0.02， 在 测 光 红 移 误差 模型 添加 了 15 x 3=45 个 自由 度 后 ， 星 系 偏 祖 因子 2， 测 
光 红 移 区 间 的 核心 部 分 的 星系 真实 红 移 分 布 的 平均 值 z, ,和 标准 差 O ore ;的 限制 情况 与 没有 灾难 性 测 
光 红 移 误差 的 情况 并 没有 实质 性 的 差别 。 图 3 比较 这 种 情况 下 与 没有 灾难 性 测 光 红 移 误差 时 BS zoei 
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和 0, ;的 限制 情况 。 重 子 声波 振荡 和 联合 重子 声波 振荡 + 弱 引 力 透 镜 对 。 的 限制 受灾 难 性 测 光 红 移 误 
差 的 影响 主要 出 现在 相应 区 间 附 近 ( 即 z<0.5 和 z>2) ， 而 中 间 没 有 灾难 性 测 光 红 移 误差 的 区 域 则 不 受 
影响 。 总 体 来 说 ， 联 合 重 子 声波 振荡 和 弱 引 力 透镜 仍然 使 星系 偏 祖 因子 得 到 比 单 独 重子 声波 振荡 更 强 
的 限制 。 男 外 ， 重 子 声波 振荡 和 联合 重子 声波 振荡 + 弱 引 力 透镜 对 zw 和 cue 的 限制 几乎 没有 变化 。 
综 上 所 述 ， 在 加 入 F,=0. 02 的 灾难 性 测 光 红 移 误差 的 情况 下 ， 重 子 声 波 振荡 和 弱 引 力 透镜 仍然 具有 
很 强 的 互补 效应 。 


Differences between with and without catastrophic photo-z errors (Ft=0.02, f;=0.073) 


9 (Zeore)*1000/(1+Zpn) 
O (Fcore)*1000/(1+Zpn) 


图 3 比较 存在 和 不 存在 局 域 性 灾难 性 测 光 红 移 的 两 种 模型 下 重子 声波 振荡 ， 弱 引力 透镜 及 联合 重子 声波 振荡 + 弱 引 
力 透 镜 的 校正 能 力 。3 张 子 图 的 图 例 相 同 并 画 在 第 3 张 子 图 上 。 第 1 张 子 图 是 重子 声波 振荡 和 联合 重子 声波 振 
荡 + 弱 引力 透镜 对 星系 偏 祷 因子 的 限制 情况 ， 并 分 别 对 应 两 条 曲线 ， 略 高 者 为 存在 灾难 性 测 光 红 移 的 情况 (后 面 
两 子 图 同 理 ) 。 每 个 参数 的 先 验 误差 为 15%。 第 2 张 子 图 是 弱 引 力 透 镜 、 重 子 声波 振荡 和 联合 重子 声波 振荡 + 弱 
引力 透镜 对 zwei 的 限制 情况 。 先 验 误差 为 Az,,, ;=0.050,i/V1-f;。 第 3 张 子 图 是 对 0, ,的 限制 情况 。 先 验 误 
差 是 Az,, ;的 1/Y2 倍 ， 对 重子 声波 振荡 和 联合 重子 声波 振荡 + 弱 引 力 透 镜 ， 图 中 几乎 分 辩 不 出 存在 和 不 存在 灾 
难 性 测 光 红 移 的 限制 差别 。 子 图 2，3 中 弱 引 力 透 镜 的 限制 情况 在 0<z,,<0.5 和 2<zw<3 出 现 明显 的 差别 是 由 于 

先 验 误差 的 取 值 不 同 造成 的 
Fig.3 The comparison of different constrains by BAO, WL and BAO+WL between with and without localized catastrophic photo- 
z error (F,=0.02). The legends for each panel are the same and put in the 3" panel. The 1“ panel is constrains of 
galaxy biases by BAO and BAO+WL and each has two curves respectively corresponding to the one with and without 
catastrophic photo-z error respectively. For both BAO and BAO+WL, the curve which is a little higher is the one with 
catastrophic photo-z error (The same for the 2" and 3" panel). The prior for each galaxy bias parameter is 15%. The 2" panel 


is constrains of z.e ; by BAO, WL, BAO+WL. The prior is Azore ; =0.050,,./ /1-f; . The 3" panel is constrains of T oore i* 


The prior for Oooo; is Azo ;/ V2. For BAO and BAO+WL constrains, the curves with and without catastrophic photo-z errors 


core ,i core ,i 


<0. 5 and 2<z „<3 is due to 


are almost the same in the 3" panel. The obvious differences between WL constrains at 0<zm ph 


the different priors adopted 


表 3 列 出 了 加 入 灾难 性 测 光 红 移 误差 到 测 光 红 移 误差 分 布 模型 中 之 后 的 时 能 量 状态 方程 参数 的 限 
制 情况 ， 并 与 没有 灾难 性 测 光 红 移 误差 时 的 结果 做 了 比较 。 在 加 入 F, =0. 02 的 灾难 性 测 光 红 移 误差 
后 ， 总 体 来 看 ， 暗 能 量 状态 方程 参数 的 误差 并 没有 太 大 的 增加 。 特 别 是 重子 声波 振荡 的 限制 结果 增加 
量 少 于 1% 。 弱 引力 透镜 对 w M w, 的 限制 误差 分 别 增加 了 14% 左 右 ( UL+BR) 和 6% 左右 (UL+BR)。 
而 联合 重子 声波 振 功 + 弱 引力 透镜 对 w。 Mw, 的 限制 误差 都 只 增加 了 5% 左右 (UL+BR) 。 


3 总 结 及 讨论 
本 文 建构 了 在 z-z, 面 左上 角 和 右 下 角 局 域 性 分 布 的 灾难 性 测 光 红 移 误差 分 布 模 型 。 计 算 了 由 于 


不 知道 灾难 性 测 光 红 移 误差 的 存在 ( 环 =0.02) 导 致 的 重子 声波 振荡 、 弱 引力 透镜 以 及 联合 重子 声波 振 
荡 + 弱 引力 透镜 限制 瞳 能 量 状态 方程 参数 时 的 系统 偏差 ， 并 分 别 考察 左上 角 和 右 下 角 部 分 的 测 光 红 移 
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表 3 加 入 灾难 性 测 光 红 移 (下 ,=0.02) 到 拟 合 模型 中 后 暗 能 量 状态 方程 参数 的 限制 结果 
Table 3 Constrains of dark energy EOS parameters after adding the catastrophic 
photo-z error ( F,=0. 02) to the fitting model 


add cat’ R”(10°) 

UL BR BOTH UL BR BOTH 
BAO Wo 0. 11960 0. 1189 0. 11963 0. 67 0. 08 0. 70 
w, 0. 4217 0. 4193 0. 4218 0. 62 0.05 0. 64 

WL Wo 0. 097 0. 095 0. 10 10 8.0 14 

w, 0. 355 0. 345 0. 360 4.4 1.5 5.9 

Wo 0. 065 0. 066 0. 067 1.6 3.1 4.7 

BAO+WL 
w, 0. 2120 0.2119 0. 2159 1.9 1.8 5.5 


result of “add cat” 


result of “no cat” 


a. “cat” means “catastrophic photo-z error” ; b. R = 


误差 造成 的 影响 。 结 果 发 现 , 左上 角 和 右 下 角 造 成 的 偏差 有 些 是 正 负 反 号 的 (重子 声波 振荡 对 w, 35 
引力 透镜 对 w ， 联 合 重子 声波 振荡 + 弱 引 力 透 镜 对 wo) ， 特 别 在 对 弱 引 力 透镜 限制 ww 的 影响 上 ， 左 
上 角 和 右 下 角 造 成 的 偏差 大 小 几乎 相等 ， 但 符号 却 相 反 ， 从 而 导致 总 体 的 影响 反而 变 得 很 小 。 重 子 声 
波 振 荡 受 灾难 性 测 光 红 移 误差 的 影响 最 小 ， 对 于 f=0.02， 对 重子 声波 振荡 最 大 的 影响 是 左上 角 + 右 
下 角 两 部 分 对 w 造成 的 偏差 (31%)。 最 后 把 F, =0. 02 的 灾难 性 测 光 红 移 误差 加 入 测 光 红 移 误差 模型 
中 ， 计 算 增加 45 个 自由 度 后 对 星系 偏 祖 因 子 、 测 光 红 移 误差 参数 以 及 暗 能 量 状态 方程 参数 的 限制 。 
结果 在 加 入 F,= 0. 02 的 灾难 性 测 光 红 移 误差 的 情况 下 ， 重 子 声波 振荡 和 弱 引 力 透 镜 仍 然 具 有 很 强 的 
互补 效应 。 重 子 声 波 振荡 和 联合 重子 声波 振荡 + 弱 引力 透镜 对 星系 偏 祖 因 子 的 限制 受灾 难 性 测 光 红 移 
误差 的 影响 主要 出 现在 相应 有 灾难 性 测 光 红 移 误差 的 区 间 附 近 ( 即 z<0.5 和 z>2) ， 而 中 间 没 有 灾难 性 
测 光 红 移 误差 的 区 域 则 几乎 不 受 影 响 。 总 体 来 说 ， 联 合 重子 声波 振荡 和 弱 引 力 透 镜 仍然 使 星系 偏 祖 因 
子 得 到 比 单独 重子 声波 振荡 更 强 的 限制 。 另 外 ， 重 子 声波 振荡 和 联合 重子 声波 振荡 + 弱 引 力 透 镜 对 
zoe i 和 0, ;的 限制 几乎 没有 变化 。 至 于 对 暗 能 量 状态 方程 参数 限制 的 影响 ， 在 把 f=0. 02 的 灾难 性 测 
光 红 移 误 差 加 入 拟 合 模 型 后 ， 暗 能 量 状态 方程 参数 的 误差 并 没有 太 大 的 增加 。 特 别 是 重子 声波 振荡 的 限 
制 结果 增加 量 少 于 1%。 弱 引力 透镜 对 wy 和 w, 的 限制 误差 分 别 增加 了 14% 左 右 (UL+BR) 和 6% 左 右 (UL 
+BR) ， 而 联合 重子 声波 振荡 + 弱 引 力 透 镜 对 w, Mw, 的 限制 误差 都 只 增加 了 5% 左 右 (UL+BR) 。 
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Abstract: The impact of the catastrophic photo-z errors on the constraints of dark energy equation of state 
(EOS) parameters and the complementarity of baryon acoustic oscillations (BAO) and weak lensing (WL) is 
studied in this paper. Targeting LSST ( Large Synoptic Survey Telescope ) -like survey projects, we construct a 
localized catastrophic photo-z error distribution in the upper-left (UL) and bottom-right (BR) of the z-z, 


plane. For the two parts of the catastrophic distribution, we use Fisher matrix to forecast the impact on dark 
energy EOS parameter constraints by the methods of BAO, WL and joint BAO and WL (BAO+WL). If the 
existing catastrophic photo-z errors are not included in the fitting model, systematic biases on EOS parameters 
could be presented. The constraints from UL and BR do not always have the same sign. BAO suffers the least 
from catastrophic photo-z errors. For total fraction of catastrophic outliers F, =0. 02, the biggest bias is on wo for 
BAO (UL+BR) which is about 30% statistical error. But for WL or BAO+WL, biases are always several times 
the statistical errors, which shows that the F, =0.02 catastrophic photo-z errors can not be ignored. Moreover, 
the impacts of UL and BR on wo for WL are almost the same magnitude but different signs, which reduces the 
total impact. With the F, =0. 02 catastrophic photo-z errors included in the fitting model, the complementarity of 
BAO and WL is still very strong even though 45 degree of freedom are added into the photo-z error distribution 
model. Under this condition, the errors of EOS parameters do not increase much. In particular, the errors of EOS 
parameter from BAO increase by less than one percent compared to that without catastrophic photo-z errors. The 
errors of wọ and w, from WL increase by approximately 14 percent (UL+BR) and 6 percent (UL+BR), 
respectively. For BAO+WL, the errors of wu and w, both increase about 5 percent (UL+BR). 
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